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Lamellar perovskite-based copper oxide compounds display three remarkable prop- 
erties: rf-wave superconductivity, with unprecedentedly high transition temperatures (/), 
a nontrivial ('Mott') insulating state (2) and a 'pseudogap' regime of suppressed den- 
sity of states (3). P. W. Anderson (2) argued that these have a common origin as strong- 
correlation phenomena that could be understood in terms of the two-dimensional Hub- 
bard model, but the interplay between the pseudogap and the superconductivity has re- 
mained controversial, with some authors suggesting the pseudogap is a signature of un- 
usual superconducting fluctuations (4, 5) and others suggesting it is a competing phase or 
regime (3, 6). Here we use new numerical techniques (7, 8) to solve the Hubbard model 
in the dynamical cluster approximation (9), which becomes exact as cluster size tends to 
infinity. Our methods enable the study of interactions strong enough and temperatures 
low enough to that the properties of the superconducting state and its relation to the Mott 
insulator and the pseudogap can be determined for clusters large enough to be representa- 
tive of the thermodynamic limit and show that in the Hubbard model the two phenomena, 
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while both linked to proximity to the Mott insulating state, are competing phenomena. 
Superconductivity is found to occur in a dome separated from the Mott insulating state by 
the pseudogap. The emergence of superconductivity from the pseudogap regime leads to 
a decrease in the gap (i.e. the creation of new states within the pseudogap), in agreement 
with recent angle-resolved photoemission data (10). 

Determining the properties of the superconducting state of the two dimensional Hubbard 
model, defined as 



has been a central issue in condensed matter theory. The existence of a superconducting state 
was suggested by numerical and analytic studies of superconducting correlations in related 
one-dimensional "ladder" systems (11). This idea was supported by analytic calculations, orig- 
inally perturbatively about suitably defined weak coupling limits (12) and later improved by 
renormalization group methods (13). Early attempts to access the superconducting phase nu- 
merically led to inconsistent results (14-17), but the development of cluster dynamical mean 
field (DMFT) methods, which become exact as cluster size iV c — > oo, (9) provided a new ap- 
proach. In pioneering papers Lichtenstein and Katsnelson (18) and Maier et al. (19) showed 
that these methods lead to <i-wave superconductivity in four-site cluster approximations. Fur- 
ther studies of the superconducting state of four site clusters were performed (20-23), but we 
now know (9, 20, 24) that the four-site cluster is in many respects an outlier, both because the 
isolated four site cluster (but not larger clusters) has an RVB singlet ground state which dom- 
inates the physical properties and because the four site cluster lacks the momentum resolution 
needed to exhibit the sector selective transition which is the expression of the pseudogap in 
this model (24). Landmark studies of the pairing susceptibility of larger clusters performed by 
Maier and collaborators (25) demonstrated that d x 2_ y 2 superconductivity is a property of the 
infinite cluster size limit. However, up to now the superconducting properties of these larger 
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clusters have been studied only at weak coupling (where the pseudogap phase was not present) 
and via analysis of the normal state pairing susceptibility (which precludes analysis of the su- 
perconducting state) (25, 26) so the question of the interplay of the pseudogap, the Mott phase 
and the superconductivity has remained open. 

Our studies use the dynamical cluster approximation (DCA) version of dynamical mean 
field theory (9). We use electron dispersion e p = —2t [cos(p x ) + cos(p y )] — At' cos(p x ) cos(p y ) 
(with the t permitting a breaking of particle-hole symmetry). We study a range of interactions 
< U < 7t and temperatures as low as t/60, enabling access to the low-T properties of the 
superconducting state and to interactions strong enough to drive a Mott transition. Submatrix 
update techniques (7) are crucial; also, in some parameter regimes care is required in converging 
to the superconducting state (see supplementary material). We present results for clusters with 
N c = 4, 8, 16 sites, large enough to distinguish generic behavior from that specific to particular 
clusters (24). 

We study the paramagnetic phase of the model, allowing for all possible spin singlet super- 
conducting order parameters. Consistently across all clusters and in agreement with previous 
results (12, 13, 25) we found that the Hubbard model possesses a d-wave superconducting state, 
with a typical transition temperature ~ t/40 « 100/<" (using a t m 0.3eV representative of the 
Cu02 superconductors). At the temperatures accessible to us, superconductivity occurs only for 
interactions larger than a critical value and dopings sufficiently close (< 20%) to half filling, 
but persists up to the largest interactions (U/t — 8) we have investigated. 

The right panel of Fig. [Qpresents the phase diagram determined from a comprehensive sur- 
vey of parameter space for N c = 8, which previous work (24) shows adequately represents the 
normal state physics of the model. Studies of selected U and doping values in the computation- 
ally much more expensive N c = 16 site cluster confirm that the physics found for N c = 8 is 
generic. The phase diagram in the absence of superconductivity has been determined (24). One 
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key feature is the presence of a fully gapped ("Mott") insulating phase for density n = 1 and 
U > U c ~ 6t indicated by the vertical heavy green line. The Mott phase is separated from the 
normal Fermi-liquid phase by a pseudogap regime characterized by a suppression of density of 
states near the (0, it), (tt, 0) regions of the Brillouin zone, with states near the zone diagonal re- 
maining unaffected. The onset of the normal state pseudogap is indicated by the dashed purple 
line in Fig. [T] 

The red shaded region in Fig. Q] indicates the parameters where superconductivity is found: 
at the temperatures accessible to us it 'turns on' as the pseudogap regime is approached, is 
strongest at the normal state pseudogap boundary, and is suppressed as parameters are tuned 
deeper into the normal state pseudogap regime. The superconducting regime is separated from 
the Mott insulating state by a non- superconducting pseudogapped regime (blue shaded region in 
Fig. [T]). The same physics is found if interaction strength is varied at doping n = 1. The lower 
left panel shows the strength of the superconducting order parameter as a function of interac- 
tion strength. Although the critical interaction strength for superconductivity onset varies with 
cluster size, in the larger clusters the superconductivity is tied to the vicinity of the pseudogap. 

The situation is somewhat different in the four site cluster, which lacks the momentum 
resolution to exhibit the momentum-space dependent density of states suppression which is the 
essence of the pseudogap phenomenon (24). In this cluster superconductivity extends all the 
way to the boundary of the Mott phase as shown in the upper left panel of Fig. [Q 

Fig. |2] presents the superconducting transition temperature and gap values estimated as de- 
scribed in the supplementary material. The transition temperature also has a dome-like behavior, 
with the highest transition temperature occurring near the onset of the normal state pseudogap 
(insets of Fig. ©, whereas the gap monotonically increases from high to low doping or low to 
high interaction. We find 2A/T C ~ 7.5 — 8 in the region outside the pseudogap (about twice 
as large as the BCS value 2A/T C ~ 3.5), and becoming rapidly larger within the pseudogap 
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regime as the endpoint of the superconducting regime is approached. Note in particular that the 
superconducting transition temperature and anomalous expectation value (see inset to Fig. [Q) 
vanish for a range of U < U Mott even though the gap continues to increase up to the endpoint 
of the superconducting phase. 

Fig. [3] presents the frequency and temperature dependence of the density of states. The up- 
per panel shows spectra representative of dopings higher than, or interactions weaker than, the 
values which maximize T c , so that superconductivity emerges from a relatively conventional 
normal state. The spectra are consistent with expectations from standard theory (29): the on- 
set of superconductivity is associated with a suppression of density of states at low frequency 
and with the formation of density of states ("coherence") peaks at slightly higher energies. We 
define the superconducting gap as half of the peak to peak distance. The behavior is conven- 
tional in the sense that the area in the coherence peaks comes mainly from the states removed 
at |a>| < A (see further discussion in the supplementary material). The gap amplitude develops 
very rapidly with temperature: only at the temperature closest to T c is the peak to peak splitting 
appreciably different from its value at the lowest T. 

The situation is quite different when superconductivity emerges from the pseudogap regime. 
Representative spectra are shown in the lower panel of Fig. [3] The normal state pseudogap is 
visible at T > T c as a suppression of the density of states at low frequencies with a broad gap 
structure at higher frequencies. The normal state pseudogap maxima are marked by dashed 
lines. The development of superconductivity is characterized by the formation of coherence 
peaks at energies below the pseudogap, i.e. by a decrease in gap magnitude as the supercon- 
ducting state is entered. Furthermore, most (typically more than 50% ) of the spectral weight in 
the coherence peak is drawn from frequencies greater than than the superconducting gap ener- 
gies. Finally, as shown in Fig. [2l the superconducting gap in this regime remains large although 
the transition temperature is driven to 0. This behavior is consistent with recent experimental 
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reports (10) that in underdoped cuprates the emergence of superconductivity out of the pseudo- 
gap regime is associated with the formation of new states at energies lower than the pseudogap 
energy and that the superconducting gap is tied to the pseudogap. 

Insight into the frequency scales relevant to the pairing interaction comes from the imaginary 
part of the real-axis anomalous self-energy shown in Fig.|4] Within our numerical precision the 
only important structure is a sharp peak at the very low frequency Cl p ~ 0.15£. This scale 
is found for all of the parameters we have examined, and is consistent with a spin-fluctuation 
origin for the superconductivity (26). We have not found evidence for the contribution from 
higher frequency (u ~ U) scales reported by Ref. (22) (see also (30)). The difference may 
arise from the different interaction values studied, from the precision of our real frequency 
spectra or from the special properties of the four-site cluster, for which the frequency scale is 
clearly larger (see supplementary material). 

The existence of d x 2_ y 2 superconductivity in the Hubbard model has been previously estab- 
lished. We find, robustly over a range of cluster sizes, interaction strengths, and carrier concen- 
trations, that the superconducting phase occurs within a superconducting dome bounded on one 
side by proximity to the Mott insulating phase and on the other by the weak coupling/large dop- 
ing regime. Superconductivity is characterized by a low pairing energy scale, is maximized in 
the vicinity of the pseudogap onset but competes with the pseudogap and is driven to zero deep 
within the pseudogap regime. When superconductivity emerges from the pseudogap regime it 
fundamentally reconstructs the density of states, leading in particular to new states within the 
pseudogap. 

Our methods are restricted (at least in the present state of development) to models with lo- 
cal density-density interactions and to situations where the fermionic sign problem, which gets 
exponentially worse with increasing cluster size and decreasing T, is not too severe. In the 
context of the two-dimensional Hubbard model the practical restriction is to N c < 16 and inter- 
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action U < 7. These values are large enough to enable access to the doped Mott phase. Within 
these constraints, investigations of the effect of second neighbor coupling, of the frequency de- 
pendence of the anomalous self energy, and of the two particle (e.g. Raman) spectra should 
be feasible. Also, a significant difference between our calculations and experiment is that we 
find a larger anomalous Green function on the electron doped side. Inclusion of long-ranged 
antiferromagnetism and also extension of our results to the 'three-band' copper oxide models is 
needed to understand these issues further. 
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Figure 1 : Superconducting phase diagram of the two-dimensional Hubbard model in plane 
of interaction strength U and carrier concentration x computed using the 8-site (right panel) and 
4-site (left upper panel) DCA dynamical mean field approximation at temperature T = i/40 
with t'/t = 0. Dashed line: location of the normal state pseudogap onset. Circles (red and light 
green) and red shading indicates the superconducting region, squares (black) and no shading the 
non- superconducting Fermi liquid, diamonds (blue) and blue shading the nonsuperconducting 
pseudogap region (eight-site only, right panel) and triangles and heavy solid line (dark green) the 
Mott insulating region at n = 1 and U > U c . Light green circles denote the points analyzed in 
Fig. [3] Lower left panel: Anomalous equal-time correlator (cx=( 7r ,o)t c .RT=(-'7r,ou) as a function of 
interaction for clusters of size N c = 4, 8, 16, as a function of U /U pg . Dashed line and pseudogap 
scale U pg = 4.2 (4-site, circles), 5.6 (8-site, squares), and 3 (16-site, diamonds) determined as 
in (24, 27, 28). 
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Figure 2: Gap to transition temperature ratio 2A/T C (main panel) computed using the 8- 
site DCA approximation both by varying U at x — (open symbols, lower axis, red color, 
Umou — 6.4t) and by varying x for U = 6t (filled symbols, upper axis, black). Left inset, 
squares: doping and interaction dependence of transition temperature for same parameters, 
showing superconducting dome. Right inset, diamonds: doping and interaction dependence of 
gap 2A/t. Arrows: onset of normal state pseudogap. 
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Figure 3 : Analytically continued spectral function computed at U — 6t for antinodal sector 
showing temperature evolution of gap structure for a typical optimally doped / overdoped state 
(x = 0.076, upper panel) and underdoped pseudogap state (x = 0.034, lower panel). Solid 
lines: superconducting spectral function. Dashed lines: normal state spectral function, obtained 
for T = t/60 by artificially suppressing superconductivity. Arrows: superconducting gap size 
A. Dotted lines: pseudogap size at T = t/30 (obtained from maximum in spectral function). 
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Figure 4: Imaginary part of self energy calculated for particle-hole symmetric Hubbard model 
atn = 1,U = 5.8, andT = t/60 obtained by directly continuing the Matsubara axis self energy. 
Normal component: black solid line. Anomalous component: dashed red line. 
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1 Numerical procedure 



1.1 Method 

We study the two dimensional Hubbard model (Eq. 1 of main text) using the dynamical cluster 
approximation (DCA) version of cluster dynamical mean field theory (i), with our recent imple- 
mentation of a Nambu (superconducting-state) version of the numerically exact continuous-time 
auxiliary field (14, 24) quantum impurity solver with submatrix updates (25) based on the open 
source ALPS (2-4) libraries. 

To describe the DMFT self-consistency process it is convenient to use the Nambu matrix 
notation introducing Pauli matrices To,i,2,3 in particle-hole space so that the Green function G, 
self energy £ and mean field function Q l have normal (N) and anomalous (A) components 
(e.g. £ = TbEjv + ti£a). The mean field function is determined from 



Q Q 1 (K, iuj n ) = H{K, iuj n ) + ^ 



,tu n ) 



.k€K 



(1) 



with 



Gi att (k, iu n ) = [iu n T + (fi - e k )r 3 - Y,(K, iu n )] 1 (2) 

Here K labels momentum sectors and S is obtained from a solution of the corresponding im- 
purity model. 
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The extension to superconductivity means that all matrices are twice as large as in a normal 
state computation at the same temperature.The doubling of matrix size means that submatrix 
update techniques (5) are crucial for accessing the large interaction strengths and low tempera- 
tures needed in this study. Critical slowing down means that typically more than 50 iterations 
are needed to converge near the superconducting phase transition. 

1.2 Convergence procedure 

Particularly in the pseudogap regime, the onset of superconductivity drastically changes the low 
energy electronic structure. A consequence is that adding a superconducting component to a 
converged non- superconducting solutions leads to very slow convergence, and care is required 
in converging to the superconducting state. We have found that the most stable procedure is to 
begin at a relatively high temperature (e.g. (3 = 10 /t) and introduce a pairing field r]i{k) = r\\^k 
via the replacement G(k, iuj n \ r/i) = [iuj n T + (fi — e fc )r 3 + r]i(k)ri — S] -1 , with e.g. 4> k = 
cos k x — cos k y for d-waxe superconductivity, and r]i typically O.lt. Retaining the pairing field we 
obtain converged solutions G(K, iu n ; r/i) first at the initial temperature, then, using the solution 
at the initial temperature as a seed, at the desired range of lower temperatures. We remark that 
the sign problem for large r]i is much less severe than at rj — 0, so these computations are not 
inordinately expensive. 

Then, at each temperature, using the converged G(K, iu n ; r]i) as a seed, we set 7/1 = 
in the self-consistency condition and continue iterating until convergence is reached. At se- 
lected points we check the solution by taking the putatively converged self energy, dividing the 
anomalous part by a large number (typically 20), and verifying that under further iterations the 
solution converges back to the one previously found. 

2 Transition temperatures and gap magnitudes 

We estimated the transition temperature by performing a linear extrapolation of the square of 
the equal-time anomalous Green's function G^u^, using the three highest temperatures at 
which Ga > 0.1. We cross-checked the transition temperature estimates by computing the 
temperature at which the inverse of the normal state pairing susceptibility vanishes. 

We obtained the superconducting gap by analytically continuing the Gn(t) and reading off 
the distance between the quasiparticle peaks. We cross-checked the gap values by analytically 
continuing the self-energies (see below). 

The crucial finding, namely that the gap is smaller in the superconducting than in the pseu- 
dogap state may be seen directly from our imaginary time data. Eq. @] implies that 



In a system with a gap A Eq. [3]implies G(r = (3/2) ~ e A / 2T . Fig. \T\ shows that at (3 = 60 and 
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U = 5.8£ the G(t) computed with superconductivity suppressed indeed lies lower than that in 
the superconducting state, indicating that the normal state has a larger gap than superconducting 
state. 




Figure 1: Imaginary time Green function (normal, local component) measured in quantum 
Monte Carlo procedure for normal states (dashed lines) and superconducting states (solid lines) 
for n = 1 and U = 5.0t ( upper traces, red online) and U = 5.8t (lower traces, blue online). 



3 Analytic Continuation 

3.1 General observations 

We perform quantum Monte Carlo calculations, obtaining Green functions and self energies as 
functions of imaginary time. Densities of states and excitation spectra are obtained by inverting 
the relation between the imaginary-time Green's function G or self energy X and the associated 
spectral function A (for G) or S (for S) which is 

G(r) = / du— — ,-A(u;). (4) 



3 



The inversion of Eq.|4]is an ill-posed problem, because the kernel in Eq.@]has many very small 
eigenvalues, so its inverse has many large ones, implying that small (statistical) fluctuations in 
the input data G cause large fluctuations in A. To find a solution we employ the maximum 
entropy continuation method (<5). 

Accurate knowledge of errors in the input data is crucial for a reliable continuation. We 
estimated these errors from eight consecutive iterations of the converged solution. We assumed 
the covariance matrix to be diagonal in frequency, based on earlier work in the normal state (7). 
We assessed the quality of our continuations by verifying that the back-continuation Eq. [4] did 
not contain systematic deviations from G within the errors bars of the data. Further confidence 
comes from the fact that the results do not change significantly as the precision is increased 
and that the change in the spectral functions is gradual and systematic as parameters (U, T, (n)) 
are varied. However, uncertainties remain in the continued quantities. Data near to = are 
in general more reliable and reproducible than data at high frequencies, and our estimates of 
energy gaps and of densities of states at frequencies up to a; ~ 0.5t are robust to choices of 
model function and input data, and show reasonable trends between calculations. However, 
the continuations in particular contain small amplitude long-period oscillations which account 
for part of the differences seen between curves in Fig. 3 of the main text at a; > 0.5t and 
which complicate for example the problem of obtaining reasonable estimates for the difference 
between two spectral functions needed for the spectral weight integral / (Eq.(6]). 

3.2 Continuation of self energies, superconducting state 

Our previous experience (8) is that in situations where a gap is present, continuation of the spec- 
tral function leads to unphysical broadening of the gap edge and that more accurate estimates 
can be obtained by continuing the self energy and then reconstructing the Green function from 
the Dyson equation. However, self energy continuation in the superconducting state is subject 
to complications. 

Although analytical continuation of the anomalous self energy proceeds from Eq. |4j T, A is 
an odd function of frequency, so the issue of normalization must be handled differently. We 
rewrite the Kramers-Kronig relation in terms of S = Y!' a (uj)/uj as 



As Eq. [5] indicates, S is normalized to the zero frequency value of the anomalous self energy, 
which we obtain by a quadratic extrapolation of l/ZU(«w n ) to cu n = 0. A very accurate ex- 
trapolation is needed because an incorrect normalization will lead to spurious features in the 
continued function, either very near zero or at high frequency. 

The results of the continuation are shown in Fig. [2] and Fig 4. in the main paper. The 
upper panel of shows the imaginary parts of the normal and anomalous self energies obtained 
by separately continuing the Matsubara axis self energies obtained at n — 1 and 5.81 One sees 
that both the normal the self energies are dominated by poles at a low energy Qp' A . The same 
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phenomena can be seen directly from the Matsubara axis data, without continuation. Fig. [3] 
shows the Matsubara axis frequency dependence of the anomalous term in the electron self 
energy. We see that the behavior is very close to the ~ 1/ (u^ + coq). 

In these continuations the pole strengths are comparable (integration, not shown, over the 
frequency range 0.1 < cu < 0.2 gives 0.32t 2 for both curves) but we see that the pole posi- 
tions are not quite identical f2p ^ Vtp and the widths are substantially different with the pole 
in the anomalous self energy being narrower. However, not all of the differences between the 
continued self energies are physical. In particular, the positivity of the density of states (imag- 
inary part of trace of Nambu Green function) implies a delicate relation between the normal 
and anomalous self energies which may be violated by the inevitable errors of analytical con- 
tinuation. In fact the density of states reconstructed from the self energies shown in Fig. [2] has 
an unphysical sign change in the vicinity of the pole frequency, arising mathematically from 
the slight difference in pole position combined with the relative narrowness of the peak in E A . 
This problem does not occur in standard phonon-mediated superconductors because the pole 
strength in anomalous part of the self energy is small and the broadening is large. 
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Figure 2: Imaginary part of the normal (r , black solid line) and anomalous (ti, red dashed line) 
components of the self energy calculated for particle-hole symmetric Hubbard model at n — 1, 
U = 5.8, and T = t/60. Upper panel: self energies obtained by directly continuing Matsubara 
axis self energy. Lower panel: self energies obtained by continuing S ± = E^ ± E^ 
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We have not found a reliable procedure for ensuring that analytical continuation of the 
self energies respects the relation between ~E N and implied by the positivity of the density 
of states, but in the particle-hole symmetric case (half filling, t — 0) progress is possible 
because the Nambu Green function is diagonalized at all frequencies by the combinations G± = 
Gn ± Ga- In this case we continue £± = ± and reconstruct the normal and anomalous 
components from the sum and difference. Results obtained for n = 1 and U = 5.8t are shown 
in the lower panels of Fig. |2] At frequencies cu > 0.2, the inevitable errors in the continuations 
(in particular the small amplitude long period oscillations mentioned in the previous section) 
lead to a small unphysical (negative sign) component in S A but the basic structure (a large pole 
at oj ~ 0.15t) survives and the small negative values have a minimal effect on the resulting 
densities of states. 




n 



Figure 3: Main panel: frequency dependence of anomalous term in electron self energy com- 
puted for the 8-site DCA approximation at T = t/40, U = 7t, and densities indicated. Left 
inset: comparison of frequency dependence of anomalous self energy for 4 and 8-site cluster 
approximations computed at U = 6t,T = t/40 and x = 0.08. Right inset: product c^E^io^) 
demonstrates decay of anomalous self energy for U = 5t and x = 0. 
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4 Transfer of spectral weight between superconducting and 
normal states 




Energy [t] Energy [t] 



Figure 4: Heavy solid line (green online, left vertical axis): spectral weight transfer integral 
Eq. [6] computed as described in the text for U = 5.0t (left panel) and U = 5.8t (right panel). 
Light solid and dashed lines: superconducting and normal densities of states, divided by 15 so 
as to appear on the same scale. 

Inspection of Fig. (3) of the main text suggests that when superconductivity emerges from 
a Fermi liquid state, the spectral weight in the superconducting coherence peak is transferred 
up from energies below the superconducting gap, whereas when superconductivity emerges 
from the pseudogap state the weight in the superconducting coherence peak is transferred down 
from higher energies. Making this statement more precise is difficult because in the general 
(non-particle-hole symmetric) case the magnitude of the weight transfer is close to the limits 
imposed by the uncertainties in our analytical continuations. However, the higher- precision 
data obtainable in the particle-hole symmetric case (which does not have a sign problem) along 
with the sharper spectra obtained by our ability to continue £± in this case permits progress. 
We construct the real-axis Ejv.A from £± for n = 1 and U = 5.0 (non- superconducting state 
is Fermi liquid) and U = 5.8 (non- superconducting state has a pseudogap) as described above, 
compute the Green function by inverting the Nambu-Dyson equation and then compute 

I (to) = dx (JmGs(x) - ImGjvO)) • ( 6 ) 
Jo 

The long-period oscillations discussed above lead to an I which fluctuates between ~ ±0.02 
for tu > 0.4i; values of |/| < 0.02 are thus indistinguishable from within our accuracy. 

The heavy solid line (left axis, green online) in the two panels of Figs. |4] shows the re- 
sulting I. In both cases, the difference is initially negative, reflecting the sharp gap in the 
superconducting state. At the gap onset the difference rises sharply, reflecting the weight in the 
superconducting coherence peak. In the left panel (overdoped case) the weight difference tends 
rapidly to zero as energy is raised above the gap energy, reflecting the fact that the weight in the 
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superconducting coherence peak is drawn mainly from states pushed up from low energy by the 
formation of the gap. In the right panel (underdoped case), the difference becomes strongly pos- 
itive for a range of cu > A, illustrating the fact that in this case the superconducting coherence 
peak is formed in large part by pulling weight down from higher energy states. One sees that 
at higher frequencies / oscillates about zero. The difficulties discussed above set an absolute 
limit of about ±0.02 on our ability to determine the difference: values smaller than this are not 
distinguishable from zero within our accuracy, so we believe that for frequencies bigger than 
about 3A the difference integral is zero within our errors. 
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